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Abstract 

One  measure  of  the  controllability  of  a  linear  system  is  the  distance  from  the  pair  (A,B)  to 
the  nearest  uncontrollable  pair.  We  provide  an  easily  computable  and  nearly  attainable  lower 
bound  on  this  distance  in  terms  of  the  "staircase  form"  of  the  pair  (A,fi). 

Keywords:  Controllability 


Introduction 

Let  A  be  an  n  by  n  matrix  and  B  an  n  by  m  matrix.  The  pair  (A, 5)  is  called  controllable 
if  the  matrix  C{A,B)  ^  [B\AB\A^B\  ■  ■  ■  |A"~'fl]  has  full  rank;  otherwise  it  is  uncontroll- 
able. One  measure  of  the  controllability  of  (,A,B)  is  the  distance  from  (A,B)  to  the  nearest 
uncontrollable  pair  (A  +&A,B  +  bB):  the  distance  is 

|18A,8fi||^=  (5:|6A,;;P  +  SISB,;!')'" 

U  'J 

if  (A  +bA,B  +  bB)  is  uncontrollable  and  (A  +  hA'  ,B  +hB')  is  controllable  whenever 
||8A',8S'||f  <  ||8A,85||f.  We  will  denote  this  distance  by  dist(A,S).  dist(A,fi)  indicates  how 
much  A  and  B  may  be  perturbed  without  making  the  system  uncontrollable. 

Several  other  workers  have  attempted  to  characterize  dist(A,S).  Eising  [3]  has  shown 
that 

dist(A  ,fi )  s  inf  (T^,„  [B  \A  -  XI]  (1) 

where  (x^i^(X)  is  the  smallest  singular  value  of  X.  This  result  shows  immediately  that 
dist(2A2*,C5)=  dist(A,S)  for  any  unitary  matrix  Q.  Boley  and  Lu  [2]  have  investigated  the 
"staircase  form"  of  the  system  (A,fl),  i.e.  a  transformation  to  an  equivalent  system 
(QAQ*,QB)  where  C  is  a  unitary  matrix,  and  in  case  B  has  only  a  single  column,  QB  has  only 
its  first  entry  nonzero  and  QAQ*=H  is  upper  Hessenberg.  If  B  has  more  than  one  column 
there  is  a  "block"  version  of  this  form:  at  most  the  first  m  rows  of  QB  are  nonzero  and  QAQ* 
is  a  conforming  block  upper  Hessenberg  matrix  (more  on  this  below).  Van  Dooren  [5] 
describes  an  algorithm  for  computing  this  staircase  form.  This  algorithm  uses  the  norms  of 
the  subdiagonal  blocks  ||//,  ,_i||  to  determine  whether  (A,B)  is  controllable.  Note  that  by  set- 
ting any  W,,-i   to  zero  (H,QB)   becomes  uncontrollable,  so  that  min||i/,  ,_i||   is  an  upper 

bound  on  dist(A,B).  Boley  and  Lu  [2]  have  shown  that  minHW,  ,_i||  can  severely  overestimate 

dist(A,B),  and  have  used  this  fact  to  criticize  the  staircase  algorithm  as  a  means  of  estimating 


dist(A,B). 

In  this  note  we  give  a  lower  bound  on  dist(A,B)  in  terms  of  the  product  of  the 
•''minC^i.i-i)-  W^  show  by  means  of  an  example  that  this  bound  is  nearly  attainable.  Thus  the 
staircase  algorithm  can  be  used  to  give  a  lower  bound  on  dist(A, fi)  as  well  as  an  upper 
bound. 


Main  Result 

We  prove  the  theorem  in  detail  when  B  consists  of  a  single  column;  the  multiinput  case 
is  analogous. 

Theorem  1:  Suppose  fi  is  a  single  column.    Assume  without  loss  of  generality  that  ||fl||=l, 
||A||  =  sup  ||Ax||/||x||  =  1,    and    that    iA,B)    is    in    staircase    form    as    defined    above.     Let 

n 

P  —  Yl   1^1.1-1 1  t)£  the  magnitude  of  the  product  of  the  subdiagonal  entries  of  A.  Then 

1=2 


dist(A,B) 


n(b(n)r 


where 


b(n) 


1  +  — 
n 


3  +  5^ 


+  4 


Remark:  fc(2)==2.672  and  quickl/ decreases  towards  its  limit  (3+5^)/2=2.618  as  n  increases. 

Proof:  By  choosing  a  diagonal  matrix  D  with  Is  and  -Is  on  the  diagonal,  DAD  can  have  subdi- 
agonal entries  with  arbitrary  sign,  so  we  assume  without  loss  of  generality  that  all  A,  ,_i>0. 
Similarly,  we  can  assume  that  Bi  =  l.  Eising's  result  (1)  implies  that 
dist(A,fi)  =  crn,in[B  |A  -X/]  for  some  fixed  X.  To  bound  this  quantity  from  below,  we  need 
the  following  fact:  a^^^^lX  \Y]  s  CTmin[^]-  This  follows  from  the  fact  that  the  eigenvalues  of 
[X  \Y]-[X  \Y]*  =  XX*  +  YY*  are  all  at  least  as  large  as  the  eigenvalues  of  XX*  [4].  Now  let 
[X  ly]  =  [fi|A-\/]  where  X  consists  of  B  and  the  first  n-1  columns  of  A-\I.  Thus 
dist(A,fi)  &  o'min(^)  =  ll^'ll''  ^^'^  ^^  need  to  bound  ||X~^||  from  above,  which  we  do  as 
follows. 

X  is  upper  triangular  with  Xii  =  l  and  Xjj  =  Ajj-i.  The  superdiagonal  entries  are 
bounded  by  lXjj  +  ,|  =  lAjj-X|  <  1+  |X  |  and  |X,j|  =  |A,j-i|  <  1  for  j>i  +  l.  It  is  tedious 
but  straightforward  to  show  that  if  we  construct  the  matrix  Z  where 


Zv  = 


0 

-1-|X| 

-1 


if  i  >j 
ifi=J 
if  i=y-l 
if  i<7-l 


then    Z    '    is    a    nonnegative    upper    triangular    matrix    with    Z^'  ^   l^i;'I-    This    implies 
||Z~'||  a  ||X"'||.  Since  ||Z~'||  <  n-max  \Z~j^\,  it  suffices  to  bound  the  largest  entry  in  Z~^ 

Again,  it  is  tedious  but  straightforward  to  show  that  the  largest  entry  in  Z~'  is  Zi„^,  the  entry 
in  the  upper  right  corner. 

For  convenience  let  j,  =  A,  +  i_,,  i'o=l,  and  a  =  1+  |X|.  We  can  define  the  following 
recurrence  for  the  entries  of  the  last  column  of  Z~': 


^.-il 


^n-2Yn 


Z-li.„  =  s-l2ia-Z-„^) 

"Z-Jn     =    i7-i(a-Z;  +  i„   +   Zj  +  2.n   +      •   •   •     +   Z„„) 

Letting  Yj  =  ZJ^'-(s„_i  •  •  •  j^_])  we  may  rewrite  this  recurrence  as 

Y„  =  1 

Y„-,  =  a 

Yj  =  aYj  +  ^  +  •^;^>+2  "•■  ^j^j  +  xYj+i  +    •  •  •    +  J/  •  •  • 
Since  0<j,sl,  Yj  is  bounded  above  by  Yj  where 

?„  =  1 

?„_,  =  fl 

P,-  =  aYj^,  +  P,+2  +  r,  +  3  +     •  •  •    +Y„ 
We  can  solve  this  last  recurrence  explicitly 

[l-a  +  ia-^-la+Sy^]    {  a  +  \- {a^-2a +5)"^"'' 
{      2(a^-2a+5)"2      J    i  2"    ""        J 

Since  a  =  1+  |X.  |a  1,  the  first  term  in  this  formula  dominates  the  second  so 


Z-,n'    = 


p 


p 


2_ 
P 


"a-l+(fl^-2a+5)'^^]  .  (  a  +  1+ (a^-2a  +  SV'^^ 
2(a2-2fl+5)'2      J    [  2 


n-l 


P     [  2 


or 


dist(A,B)  > 


nZin^ 


n 


[|X|  +  2+(|X|2  +  4)i'M 


fci(l^l) 


(2) 


This  bound  decreases  as  |\|  increases,  so  we  need  another  bound  that  increases  as  |\| 
increases  in  order  to  get  a  bound  independent  of  |X.|.  To  do  this  we  use  the  same  trick  as 
before:  dist(i4,B)  =  crn,i„[fi  |A -X/]  s  crn,in[A -X./].  As  before,  we  need  to  bound 
||(A-\/)~'||  from  above  to  bound  dist(A,B)  from  below.  If  |X.|>1,  then  since  ||A||=1 

||(A-\/)-i|i  =  ||x-H/-^/M-Nl  ^  \^-'\  I  (i-||A/x||)  =  1  /  (|x|  -  1) 

Thus  we  have  our  second  lower  bound  on  dist(A,B) 

dist(A,fi)  s:   |X|-1  =  fc2(l>^l)     if    1^   l>^l    • 

The  lower  bound  fci(|X|)  is  a  decreasing  function  of  |X|  and  ^aCl^l)  is  an  increasing 
function  of  |\|.  Since  i>2(l)<fri(l)  and  b2{2)>b^(2),  b^{\^)=b2{K\)  for  some  1  <  \,  <  2. 

Solving  for  Xj  we  get 

2 


bx{\x)=  ^ 


X,  +  2+(Xf  +  4) 


1/2 


=    Xi-1   =    b2{\l) 


Since        ^2(^1)         has        slope         1,         X.]         is         clearly         bounded         above 
1+Z>i(l)  :s  l+n~'(2/(3  +  V5))"  ^  Xj.    This  implies  that 


by 


dist(A,B)  s  fciCXi)  =  ■^• 


3  + 


3  +  5* 


+ 


1  + 


.3  +  5-J  J 


n\2  11/2 

+  4 


as  desired,  n 


Now  we  give  an  example  that  shows  that  dist(A,B)  can  in  fact  be  as  small  as  p .  Let  B  be 
the  first  column  of  the  n  by  n  identity  matrix  as  before  and  A  bidiagonal  with  Ajj=  —1,  j<n, 
A„„  =  0,  and  A^j-i  =  e^-cl.  Then  ||A||~1,  and  dist(A,B)  <  a^,^[B  \A]  =  a^inCX),  where  X 
consists  of  B  and  the  first  n  —  \  columns  of  A,  since  the  last  column  of  A  is  zero.  X~'  is  easily 
calculated  and  its  largest  entry  Xf„'  is  seen  to  be  H^/'"/'''  ^°  dist(A,B)  <  p  as  desired. 

j 
Also,  Boley  [1]  gives  an  example  where  dist(A,B)  decreases  as  2   "  with  constant  p,  so  the 

factor  b{n)~"  in  the  bound,  while  not  necessarily  attainable,  does  reflect  the  behavior  of  the 
problem.  If  either  the  minimizing  X  in  (1)  is  close  to  zero,  or  if  each  subdiagonal  entry  A,  +  ]  ,• 
is  very  small,  one  can  show  that  fc(n)~"  can  be  replaced  by  a  factor  close  to  2~".  Note  also 
that  when  p  is  small  Xj  in  the  proof  is  close  to  1  and  the  b{n)~"  factor  approaches  its  limit 
2.618~". 

To  state  the  general  result,  we  need  to  define  staircase  form.  Suppose 


B  = 


0 
0 


and     A  = 


An 

A21 

0 


'ij 


0  A, 


(3) 


where  A  is  block  upper  Hessenberg  as  shown,  fi]  has  as  many  rows  as  An.  Bi  has  full  row 
rank,  and  each  subdiagonal  block  A,  ,_i  has  full  row  rank.  Then  we  say  (A,B)  is  in  staircase 
form.  As  mentioned  above,  there  is  a  unitary  matrix  Q  that  transforms  any  pair  (A,fi)  into 
staircase  form  {QAQ* ,QB)  [5]. 

Theorem  2:  Let  B  be  an  «  by  m  matrix  and  A  an  n  by  n  matrix.    Assume  without  loss  of  gen- 
erality that  1|A||=  ||B||=  1  and  that  {,A,B)  is  in  staircase  form  with  j  blocks  as  in  (3).    Let 
;■ 
i(^i)II  ^min(^i,i-i)-  Let  b  {■)  be  defined  as  in  the  statement  of  Theorem  1.  Then 
i=2 

dist(A,B) 


cr- 


j-{bij)y 

Since  generically  j=  \nlm\,  this  lower  bound  is  an  increasing  function  of  m  (for  fixed 
P)- 
Conclusions 

We  have  shown  that  the  staircase  form  of  a  single  input  control  system  can  be  used  to 
provide  a  lower  as  well  as  an  upper  bound  on  the  distance  to  the  nearest  uncontrollable  sys- 
tem. This  bound  is  easily  computed  from  the  staircase  form,  and  is  shown  to  be  nearly  attain- 
able for  some  systems.  The  gap  between  the  upper  and  lower  bounds  can  be  very  large,  and 
so  the  problem  of  accurately  and  inexpensively  estimating  the  distance  to  the  nearest  uncon- 
trollable system  remains  open. 
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